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ABSTRACT 

X-ray observations of novae, nova-1 Ike variables, and dwarf novae In 
outburst have shown them to be weaker X-ray sources than expected from models 
attributing the X-ray emission to a boundary layer formed between the white 
dwarf surface and an accretion disk. A resolution of this discrepancy is 
Important for understanding the accretion process in cataclysmic variables. 
Constraints on the X-ray emission from novae can be derived from the lack of 
detections of novae In the HEAO-1 soft X-ray sky survey. An analysis of these 
constraints suggests that the discrepancy can be resolved If novae, nova-like 
variables, and dwarf novae In outburst have optically thick boundary layers 
with kT « 15-30 eV, provided the neutral hydrogen column densities of 
clrcumstellar 0 as are > i q20, 5_i q21 GR1 - 2. Tf so , we do not observe' a strong 
hard X-ray flux (E > 1 keV) from these systems because an optically thick 
boundary layer will not be a hard X-ray source. Likewise, we do not observe a 
strong soft X-ray flux (E < 0.5 keV), even though it Is produced In the 
boundary layer, because most of this flux is absorbed by clrcumstellar gas. 

The importance c; the photoelectric absorption of soft X-rays by clrcumstellar 
gas may have been underestimated in previous work. The- weak hard X-ray fluxes 
commonly observed from novae, nova-like variables, and dwarf novae In outburst 
may be produced in an accretion disk corona. 

Subject headings: stars: accretion - stars: dwarf novae - 

stars; novae - X-rays: binaries 


I. INTRODUCTION 


Accretion models for cataclysmic variables (CVs) predict the formation of 
a boundary layer at the interface between the white dwarf surface and an 
accretion disk (Lynden-Bell and Pringle 1974). The boundary layer Is expected 
to radiate an EUV/soft X-ray flux If It Is optically thick (Pringle 1977) or a 
hard X-ray flux If it Is optically thin (Pringle and Savonlje 1979; Tylenda 
1981). Approximately half the gravitational energy released by accretion 
should be released by dissipation of rotational energy In the boundary layer 
so that boundary layer luminosities greater than 10 33 ergs s -1 are expected. 
X-ray observations of CVs, however, have shown most of them to be weaker X-ray 
sources, with X-ray luminosities less than ID 32 ergs s" 1 (Cordova and Mason 
1983). A discrepancy between the observed low X-ray luminosities and the 
luminosities expected from the boundary layer has been noted, particularly for 
the nova subclass of CVs, by Cordova, Jensen, and Nugent (1981; hereafter 
CJN), Becker (1981), and Ferland et al. (1982a). It Is difficult to explain 
the absence of a boundary layer if the standard accretion model for CVs Is 
correct. A resolution of this discrepancy is therefore important for 
understanding the accretion process In CVs. 

The purpose of this paper is to extend the CJN treatment of HEAO-1 
constraints on the X-ray emission from novae to provide a better understanding 
of the boundary layers In novae. The HEAO-1 constraints are particularly 
useful for determining whether the boundary layers in novae are optically 
thick or optically thin, and for estimating the boundary layer temperatures 
and neutral hydrogen column densities to the boundary layer. 

In §11, I summarize the predictions of the boundary layer model for the 
X-ray emission from novae. A comparison with the X-ray observations Is made 
In a review of previous suggestions that a discrepancy exists between the 


observations and the theory. In §111, I extend the work of CJN to provide 
constraints on the nature of boundary layers In novae, based on the lack of 
detections of novae In the HEAO-1 soft X-ray survey. As a result, I will show 
that hard X-rays from novae are unlikely to be produced In a boundary layer. 

I will also provide estimates of temperatures and column densities for 
optically thick boundary layers In novae. In § IV, the applicability of these 
results to other subclasses of CVs, and their Implications for the nature of 
boundary layers, clrcumstellar gas In CVs, and the hard X-ray source In CVs, 
are discussed. 

II. THE NOVA PROBLEM 

a. The Standard Model 

The standard accretion model for CVs makes the following predictions. 

Approximately half of the gravitational energy released by accretion will 
be radiated from the accretion disk. Accretion disk luminosities of 6MM /2R = 
8.5 x 10 33 (M/M 0 )(M /10“ 9 M 0 yr _1 )(5 x 10 8 cm/R) ergs s” 1 are expected, where 
M is the mass of the white dwarf, M is the mass accretion rate onto the white 
dwarf, and R is the radius of the white dwarf. The accretion disk will not be 
hot enough to be an X-ray source (Shakura and Sunyaev 1973). 

The remaining half of the energy released will be radiated from a 
boundary layer between the white dwarf surface and the accretion disk. The 
luminosity of the boundary layer should be approximately half the total 
accretion luminosity. 

If the boundary layer Is optically thick. It should radiate as a 15-30 eV 
blackbody from an annulus with Inner radius R and outer radius R+h, where 
h(«R) is the semi -thickness of the inner accretion disk (Pringle 1977). Half 
of this radiation will be occulted by the white dwarf. The observable 
luminosity from this blackbody component should therefore be approximately on 


quarter of the total /Accretion luminosity. A fraction of this radiation will 
be at energies 0.1-0. 5 keV, accessible to soft X-ray detectors. The size of 
this fraction depends on the temperature of the boundary layer and the column 
density of absorbing gas. 

If the boundary layer Is optically thin, It will radiate primarily by 
thermal bremsstrahlung from a 1-30 keV plasma (Tyjenda 1981). Half of this 
radiation will be Intercepted by the white dwarf and accretion disk and 
re-radiated at longer wavelengths. The observable luminosity from the hot 
boundary layer plasma should be approximately one quarter of the total 
accretion luminosity. Most of this flux will be a hard X-ray flux at energies 
greater than 1 keV. 
b. X-Ruy Observations of Novae 

The HEAQ-1 soft X-ray survey of 102 novae did not detect X-ray emission 
In the 0.18-2.8 keV energy band from any of them, Indicating that their soft 
X-ray luminosities are less than 10 31 ( d/1 00 pc) 2 ergs s" 1 (CJN). Einstein 
observations of 14 of these novae in the 0.1-4. 5 keV band are reviewed by 
Cordova and Mason (1983). Seven of these were detected with X-ray 
luminosities of 10 3 * - 10 32 ergs s'*. The novae detected had hard X-ray 
spectra, with most of the flux at energies greater than 1 keV. The remaining 
7 novae were not detected, indicating their X-ray luminosities were less than 
10 32 ergs s" 1 . There was no evidence for a soft X-ray component from any of 
the 14 novae observed. The conclusion from these observations is that these 
novae are weak, hard X-ray sources with X-ray luminosities less than 10 32 ergs 
s’ 1 and source temperatures greater than 0.5 keV. 

It is probable that the novae examined with HEAO-1 and Einstein 
constitute only a small fraction of the entire class of novae. Older novae, 
those whose last outbursts occurred more than 100 years ago, are probably much 
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more numerous, since the recurrence time scale for nova outbursts Is believed 
to be much greater than 100 years (Bath and Shavlv 1978), If so, there will 
be many older novae at distances nearer to the sun than are the observed 
novae. These novae have not been Identified, but some of them would appear In 
sky surveys as serendipitous X-ray sources If their X-ray luminosities were 
sufficiently high. More stringent constraints on the X-ray luminosities of 
novae can be obtained from the lack of detections of older novae as 
serendipitous sources in X-ray sky surveys. 

From the lack of detections of older novae in the HEAO-1 soft X-ray 
survey, CJN determined upper limits to the X-ray luminosity of the novae as a 
class of objects. They assumed the novae are uniformly distributed In the 
solar neighborhood with a number density n. If the typically have soft 
31) eV blackbody spectra wltt a neutral hydrogen calumrt density N^ * 10 20 cm" 2 * 
CJN derived L x < 4 x 10 29 (n. 4)" 2 / 3 ergs s 1 " 1 , where is the number density 

of novae in units of 10"^ pc" 3 . If the novae typically have hard 10 keV 
thermal bremsstrahlung spectra with N^ « 10 20 cm" 2 , they derived L x < 10 30 
(n _^)“ 2//3 ergs s"*. Becker (1981) made a similar analysis based on the lack 
of detections of older novae as serendipitous sources by the Einstein Imaging 
Proportional Counter. For the same hard X-ray spectrum considered by CJN, 
Becker derived L x < 5 x 10 29 (n .^)" 2 / 3 for older novae as a class of objects. 

The classical novae, considered as a class of X-ray sources, therefore 
appear to have X-ray luminosities much lower than expected If the boundary 
layer models for the X-ray emission from CVs are correct. This "nova 
problem", the discrepancy between the observed low L x from novae and the 
higher L x predicted by theory, has also been considered by Ferland et al. 
(1982a), who note that a resolution of this discrepancy is Important for our 
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understanding of the accretion process in CVs and the nature of the nova 
outbursts. 


III. HEAO-1 CONSTRAINTS 

The CJN and Becker analyses are limited in that a uniform distribution of 
novae was assumed and only two specific source spectra were considered. A 
more realistic treatment should account for scale heights in the distribution 
of novae in the solar neighborhood. Even more critical is the need to 
consider a wider range of source temperatures and column densities, since the 
sensitivity of the X-r^y detectors often depends strongly on these parameters. 

Another limitation is that constraints on the X-ray luminosity can be 
used to test boundary layer models only if a total accretion luminosity, hence 
a mass accretion rate, is assumed. A more useful constraint would be on the 
ratio of the boundary layer luminosity to the total accretion luminosity. 

Since boundary layer models predict that one quarter of the accretion 
luminosity should emerge as a boundary layer component, this ratio should 
be ~ 0.25, Independent of the mass accretion rate. 

In thi o section, I present an extension of the CJN work on constraints on 
the X-ray emission from novae, using results from the HEAO-1 soft X-ray 
survey, in which these limitations are overcome. A more realistic 
distribution of novae is used and a wide range of temperatures and column 
densities are examined. Rather than attempt to constrain the X-ray luminosity 
for an assumed number density and mass accretion rate, I attempt to constrain 
the fraction of the total accretion luminosity which can be radiated as an 
X-ray component solely as a function of the spectrum, temperature, and column 
density of that component. As a result, I am able to determine the types of 
X-ray spectra, temperatures, and column densities that are consistent with 
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boundary layer models. I will show that the ’’nova problem" may be resolved If 
novae have optically thick boundary layers and higher column densities than 
have previously been supposed, 
a. The HEAO-1 Survey 

The HEAO A-2 low energy detectors (LEDs) have surveyed over 95% of the 
sky In the energy bands 0.18-0.44 keV, the “1/4 keV" band, and 0.44-2.8 keV, 
the "1 keV" band. As an all sky survey sensitive to the soft X-rays expected 
from an optically thick boundary layer and the harder X-rays expected from an 
optically thin boundary layer, It Is a useful survey for examining the X-ray 
emission from the class of novae In the context of boundary layer models. 

The HEAO A-2 LED source catalogue (Nugent et al. 1983: hereafter NEA) 
uses data collected In this survey by LED!. The results Include new 
detections of 54 X-ray sources In the 0.18-2.8 keV energy band, of which 32 
remain unidentified. Of these 32 sources, 22 are detected In the 1/4 keV band 
and 13 are detected In the 1 keV band, with 3 sources detected In both 
bands. These results set upper limits to the number of novae <N X > which are 
expected to be detectable by LED1 . It Is Improbable that many of the 
unidentified sources are novae, since none of the 82 identified LED sources 
are novae. A more realistic expectation Is <N X > < 1. 

Since LED! Is a proportional counter which records photon events 
(Rothschild et al. 1979), the sensitivity of the survey Is Initially measured 
as a count rate sensrdvity. The characteristic count rate sensitivity of the 
survey in each band is determined by the criteria for Inclusion of a source in 
the catalog. In the following analysis, 4.8 counts s -1 is used as the 
characteristic count rate sensitivity of the HEAO-1 soft X-ray survey in each 
band. This Sensitivity corresponds to the 6.1a significance criterion for the 
count rate intensity which was used for including a source in the catalog 
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(NEA) . A$ a result of this survey, we may expect that there are < 1 novae 
detectable by LED1 at a count rate of 4.8 counts s** 1 . This result can be 
applied to an analysis of the constraint?. $n the X-ray emission from the class 
of novae. 

b. Method of Analysis 

Consider an X-ray source characterized by a photon spectrum dN/dE photons 
s 1 * 1 keV" 1 . The luminosity of this X-ray component Is; 

L xc * / dE[E dN/dE],* 

L xc Is the t otal luminosity of the X-ray source at al l energies, and not just 
the fraction of the luminosity In a given X-ray bandpass. The luminosity of 
the source In a given energy bandpass 1 Is; 

/ dE [E dN/dE] 

1 xc /"dE [E dN/dE] xc 1 

o 

where 3 ^ Is the fraction of L xc which Is radiated In bandpass 1. Our measured 
parameter Is Cj, the received count rate In the detector In bandpass 1. We 
can determine L-j from c-j by; 

L.j = 4nd c.jOt.j , 

where d Is the distance to the source and Is a factor which converts the 
detected count rate to the energy flux which would 'be received at the detector 
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in the absence of absorption. This factor depends on the spectral parameters 
dN/dE, kT, and N H . See NEA for a discussion of the conversion factor. 

Note that NEA call this factor dj rather than c^. 

If a nova at a distance d is radiating an X-ray competent characterized 
by dM/dE, kT, and N^ which Is detected by LED! In a bandpass 1, the fraction 
of the total accretion luminosity radiated by the X-ray component Is: 
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If we can determine a maximum distance d x at which a nova Is detectable by 
LED1, we can constrain f by: 


p 1 

D 4ird„ (4.8 counts s )cb 
f < ( JL) * 1 
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Constraints on f can be directly applied to the boundary layer models. If the 
boundary layer Is the source of X-ray emission from novae, the luminosity of 
the boundary layer Is the equivalent of the luminosity of the X-ray component 
L xc . From the predictions of the boundary layer models (§11), we then expect 
f ~ 0.25. It Is convenient to express f as f < <(>6, where 


4 > = (JL) . 4ird 2 . (4.8 counts s" 1 ) . (10" 11 ergs enf 2 count" 1 ) 
fflH 
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This separates f Into the product of two functions, one of 
which, 6 * 6{aj,$.j) * 6(dN/dE, kT, N^), Is solely a function of the source 
spectral parameters. The other function, Is not at all a function of these 
parameters. Since the white dwarf radius is determined by a mass-radius 
relation (Hamada and Salpeter 1961), $ Is a function of the white dwarf mass P 
the mass accretion rate, and the maximum distance at which a nova Is 

detectable by LEDli <j> « <j> (M,M, d x ). I win demonstrate that d x Is related 

« 

to M and M In such a way that 4 > Is sharply constrained for the range of white 
dwarf masses arid mass accretion rates expected for novae. Therefore f, unlike 
L x , can be constrained as a function solely of the source spectral parameters 
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To determine d x as a function of the white dwarf mass and mass accretion 
rate, I make two assumptions. 

I assume that there have occurred 12 nova outbursts at distances less 
than 1 kpc during a recent 88 year Interval. These estimates have been 
obtained from Duerbeck (1981). Using these estimates, the number density of 
novae which have had recent outbursts can be estimated, 

I also assume that the nova outbursts are recurrent phenomena with a 
recurrence time scale t r which Is a calculable function of the white dwarf 
mass and mass accretion rate. The outburst recurrence time scale Is needed to 
calculate the number density of all novae from the number density of novae 
which have had recent outbursts. For a derivation of T r (M,M), I have relied 
on the basic model describing the nova outburst as a thermonuclear runaway in 
the hydrogen rich white dwarf envelope (Gallagher and Starrfield 1978) and on 
the work of Nariai and Nomoto (1979) on the evolution of a hydrogen envelope 
In an accreting white dwarf. For a discussion of this work, see Nomoto 
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(1982). 

The results from the calculations of Narlal and Nomoto are represented In 
Figure 1, where the outburst recurrence time scale x r Is shown for mass 
accretion rates ranging from IQ"** M 0 yr”* to 10"? M 0 yr"* and for white dwarf 
masses of 1.0, 1.3, and 0.6 solar masses. The local number density of novae 
n, which Is proportional to x r {see Appendix), Is also shown. 

If these assumptions are correct, we can determine the maximum distance 
d x at which a nova can be detected by LED1 from 

1.1 <N„> 


where d x Is In units of kpe, <N X > is the expectation number of novae 
detectable fsy x r Is In units of years, and 


1(d) . £ dr ! Z l de [r e -C(100 + r 2 + 20roose) 1/2 /3.4] a . e -5(d 2 -r 2 ) 1 / Z n 


A description of the derivation of I(d x ) Is given In the Appendix. 1(d) Is 
plotted In Figure 2. For d«200 pc, I Is proportional to d 3 ; for d>>200 pc, I 
Is proportional to d 2 . This behavior Is a result of choosing a z scale height 
of 200 pc for the novae. Values of 1(d) are tabulated by Jensen (1982). 1(d) 

is a measure of the relative number of novae located at various distances d. 
For example, I(d « 0.1 kpc) - 4.6 x 10 -4 while I ( d = 1 kpc) = 0.1518 so that 
there are 330 times as many novae within 1 kpc as there are within 100 pc. 
The absolute number of novae located within a given distance d depends also on 
N 0 , x 0 , and x r (see Appendix). 
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Since t r Is a function of the white dwarf mass and mass accretion rate, 

we can determine d x as a function of M, M, and the expected number of 

detectable novae <N X >. An estimate of <N X > < 1 Is derived from the lack of 

detection of novae in the HEAO-1 soft X-ray survey. In Figure 3, I shew the 

dependence of d x on the mass and mass accretion rate If we expect one nova to 

be near enough to be detectable by LED1 (<N X > =1). There Is not a large 

deviation from the trend d x 2 a M. Therefore, <j>, which Is proportional to 

d x 2 / M, Is relatively insensitive to the mass accretion rate. Since 6 Is not 

• 

a function of mass and mass accretion rate, the dependence of f on M and M is 
the same as the dependence of <j> on M and M. 

The weak dependence of <f» on mass and mass accretion rate Is Illustrated 
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white dwarf mass and mass accretion rate shows that the constraints on f, the 
ratio of the total accretion luminosity which can be radiated as an X-ray 
component, are relatively insensitive to the mass and mass accretion rate. 
Since the wide range of white dwarf masses and mass accretion rates considered 
here will account for the expected masses and mass accretion rates for the 
novae, the upper limits I derive for f can be assumed to be applicable to the 
novae as a class of objects, even though most novae are not directly observed 
and in fact have not been identified. As, a class of objects, novae are 
therefore characterized by 
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where 4> max Is chosen as the maximum value of <f> for a given wi?,te dwarf 
mass. For a given mass, the mass accretion rate which maximizes 4 Is 
chosen. Thus, the value of f derived for a given white dwarf mass Is always 
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an upper limit. For M « 1 solar mass, a *=4.8 x 10” 4 ; for M = 0.6 solar 

MIC* A 

mSS ’ “ 8 - 7 X 10 ‘ 4 * for M " U3 SOlar “ ss - W " 7 - 2 x 10 ’ 4 - 

With 4> max determined, the types of source spectra consistent with 
boundary layer models are constrained by the requirement 


ct 1 (dN/dE, kT, N h ) x 

♦max' 10 -U- ergs "cm" 2 count-T ’ ^ °' 25 ' 

Two types of source spectra are considered. 

1. Thermal bremsstrahlung 2 : dN/dE ~ e“ E ^ T , appropriate for the 


2 For kT less than 3 keV, radiative recombination and line transitions in an 
optically thin plasma cannot be neglected and will have an effect on dN/dE 
(Blumenthal and Tucker 1974). and are not available for recombination 
spectra so this effect cannot be accounted for precisely. The effect on aj 
and 3 ^ will not be sufficiently large to change in any Important way the 
conclusions (§11 Ic) regarding optically thin boundary layers (Jensen 1982). 


consideration of optically thin boundary layers. Gaunt factors g(E) are taken 
from Karzas and Latter (1961). 

2. Blackbody: dN/dE ~ E*V(e^T . i) f appropriate for the consideration 
of optically thick boundary layers. 

Conversion factors a i (kT,N H ) In the 1 keV and 1/4 keV bands are 
calculated by NEA for thermal bremsstrahlung spectra (their Figure 8) and for 
blackbody spectra (their Figure 10). Note that NEA call these factors d^ 


*Wfc. “ * 


15 


rather than a^. For a fliven spectrum, the band with the smaller a^/e^ is 
chosen, since it produces a tighter constraint on f. 

C. Results of the Analysis 

The constraints on f, the fraction of the total accretion luminosity 
which can be radiated as an X-ray component by novae, are the principal 
results of the analysis. The upper limits to f as a function of source 
temperature and column density for a thermal bremsstrahlung spectrum are shown 
in Figure 5. The upper limits to f for a blackbody spectrum are shown in 
Figure 6. 

These figures illustrate f for the case of a one solar mass white dwarf 
under the assumption <N X > - 1. In this case, <j> max = 4.8 x 10"^. To obtain f 
for a specific mass and mass accretion rate, one can read fj>(M,M)/4 max from 
Figure 4 and adjust f accordingly. For example, if we choose a white dwarf 
mass of 1.3 solar masses and a mass accretion rate of 10" 8 M 0 year" 1 as 
typical, we find <f> - 7.0 x 10" 4 . Therefore log^g f should be Increased by 
0.16 over the log^f shown in Figures 5 and 6. Tables of c|>(M,M) are given in 
Jensen .(1982). Errors in the derivation of f due to errors in the estimations 
of the expected number of novae detectable by LED!, <N X >, the number of recent 
nova outbursts at distances less than 1 kpc, N 0 , and the dependence of the 
outburst recurrence timescale x r on mass and mass accretion rate, may be 
deduced roughly from the relation f a {<N X >/N 0 x r ) 0,7 . If the assumption <N X > 
= 1 is in error, it is more likely to be an overestimate. Therefore, unless 
Nariani and Nomoto seriously overestimate x r (M,M), errors in the estimates of 
<M X >, N 0 and x r should tend to overestimate f so that the upper limits to f 
derived here should be valid. 

1) Constraints on f 

For a given spectrum, temperature, and column densl"”' N^, the upper limit 
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to f, the fraction of the total accretion luminosity which can be radiated as 
an X-ray component from novae, can be read from Figures 5 ud 6. For example, 
In the case of a 30 eV blackbody spectrum with « 10 20 cm" 2 , Figure 6 shows 
that log^gf < -2.5, or f < .003. Therefore, If novae are characterized as one 
solar mass white dwarfs with 30 eV blackbody components and column densities 
of 10 20 cm" 2 , less than 0.3% of the total accretion luminosity can emerge from 
that component. If more than 0.3% emerged as a 30 eV blackbody component, we 
would expect to see more than one nova with LED1. If Is typically as high 
as IQ 21 cm" 2 , the constraint is relaxed to f < 0.24. Therefore, If 25% of the 
total accretion luminosity emerged as a 30 eV blackbody flux with * 10 CA 
cm" 2 , we would expect to detect one nova with LED1. 

Since we expect f ~ 0.25 for the*boundary layer component (§11), values 
of (kT, N^j) consistent with a boundary layer component He above the 
horizontal dashed lines marking f = 0.25 In Figures 5 and 6. It can be seen 
from these figures that f is less than 0.25 except for extreme temperatures 
and large column densities. These results can be applied to constrain the 
nature of the boundary layers in novae. 

11) Optically Thin Boundary Layers 

From Figure 5, It can be seen that if novae have optically thin boundary 
layers, very low temperatures and high column densities are required for 
consistency with the HEA0-1 constraints. For a column density < 10 22 cm" 2 , kT 
< 250 eV is required. Boundary layer temperatures this low would require low 
viscosities and very low mass accretion rates (Tylenda 1981) and would not be 
hard X-ray sources even If they could be produced. If novae have optically 
thin boundary layers hot enough to produce the hard X-ray fluxes observed from 
the selected novae in the Einstein IPC sample, column densities greater than 2 
x 10 22 cm" 2 are required to explain the lack of detections of novae by the 


HEAO-1 LED1, Unless the column densities to the nearest novae are this large, 
the HEAO-1 constraints show that only a small fraction of the rotational 
energy In the boundary layer can emerge as hard X-rays. It Is Improbable that 
the column densities are this large, since the hard X-ray spectrum would then 
show an absorption turnover at energies greater than 1 keV. Such turnovers 
are not observed In the spectra of CVs with observed hard X-ray fluxes, 
ill) Optically Thick Boundary Layers 

As noted earlier, f < .003 for the blackbody spectrum considered by CJN 
(kT = 30 eV, N^j = io 20 cnf 2 ). If novae have optically thick boundary layers, 
they must be characterized by lower temperatures and/or higher column 
densities. For a column density of 10 2 * cm -2 , boundary layer 
temperatures < 30 eV are consistent with the HEAO-1 constraints. Pringle's 
(1977) optically thick boundary layer model predicts temperatures of 15-30 
eV. Therefore, the view that novae typically have, optically thick boundary 
layers is consistent with the HEAO-1 constraints J_f the column densities to 
these boundary layers are typically > 10 20,5 -10 21 cm" 2 . 

IV. DISCUSSION 

Ferland. jit jiK (1982a) suggest that the most likely solution to the "nova 
problem" Is either a source of optical radiation in excess of the radiation 
produced by accretion, such as a reprocessed flux from a hot white dwarf, or a 
boundary layer structure greatly different from that predicted by current 
models. Sources of large optical excess can result in low ratios of X-ray to 
optical luminosity, even if the boundary layer converts a large fraction of 
the energy liberated by accretion into X-rays. Such sources have no effect on 
the ratio of X-ray to accretion luminosity, however. It is this parameter 
which is constrained by my analysis. By recasting the "nova problem" as a 
deficiency in L x /L acc , rather than L x /L opt , my analysis shows that the 
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solution cannot be a source of excess optical radiation. 

The preceding analysis of the HEAO-1 constraints on the X-ray emission 
from novae suggests an alternative solution which can be reconciled with 
current boundary layer models. The discrepancy between the low X-ray 
luminosities observed from novae and the higher luminosities expected from the 
boundary layer is resolved if novae have optically thick boundary layers with 
temperatures ~ 15-30 eV and column densities > 10 2 ^-10 21 cm" 2 . In that 
case, we do not observe a strong hard X-ray flux from novae because an 
optically thick boundary layer will not be a hard X-ray source. Likewise, we 
do not observe a strong soft X-ray flux from novae, even though it is produced 
in the boundary layer, because most of this flux is absorbed by an intervening 
column of gas, and re-radiated at longer wavelengths, 
a) The Importance of Photoelectric Absorption 

The constraints on the soft X-ray emission from novae are very sensitive 
to column densities in the range 10 20 - 10 21 cm" 2 because the cross section 
for photoelectric absorption ranges from 3 x 10" 21 cm" 2 to 10~ 20 cm" 2 for 
photons in the 1/4 keV energy band of LED!. Therefore, optical depths of 
unity occur for column densities of 1-3 x 10 20 cm" 2 . The effect of column 
density on the energy flux in the 180-280 eV band received from a 30 eV 
blackbody is shown in Figure 7. Most of the flux detectable by LED1 from a 30 
eV blackbody is received in this energy band. For a column density of 10 20 
cm" 2 , less than half of the 180-280 eV flux is absorbed, but more than 99% of 
this flux is absorbed by a column density of 10 21 cm" 2 . 

The effect of this much larger absorption is found in the conversion 
factors from received count rates to flux at the source, a. For a 30 eV 
blackbody with = 10 20 cm" 2 , a ~ 10" 11 ergs cm" 2 count" 1 , whereas for a 30 
eV blackbody with = IQ 21 cm" 2 , a ~ 10" 9 ergs cm" 2 count" 1 . A 30 eV 
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blackbody must therefore be 100 times as luminous to be detectable by LED! If 
Npj » 10 21 cm“ 2 than If » 10 2( ^ cm” 2 . For column densities greater than 10 2 * 
cm’ 2 the absorption Is so large that the probability of detecting a nova with 
LED! Is very small. For example, If novae have 30 eV optically thick boundary 
layers with as large as 3 x 10 21 cm’ 2 , the probability of detecting soft 
X-rays from a nova with LED! Is less than 1%. 

CJN did not consider column densities greater than 10 20 cm’ 2 . The 
discrepancy they found between the low soft X-ray luminosities Inferred from 
non-detections of novae and the higher luminosities expected from optically 
thick boundary layers is resolved if column densities are greater. 

b) The Circumstellar Gas 

Since the nearest novae are probably at distances < 100 pc, it is 
unlikely that interstellar column densities are greater than 10 20 cm’ 2 . If 
novae have optically thick boundary layers with kT ~ 15-30 eV, the required 
large absorption of the soft X-rays in the boundary layer flux means that 
there must be enough absorbing gas in the system to produce column 
densities > 10 20 * 5 cm" 2 . The flux from the boundary layer should have an 
effect on the ionization of this gas. Ferland _al_. (1982b) applied 
photoionization models to fit the ultraviolet and optical line spectra of the 
nova V603 Aquilae. They found better fits for models in which the line 
forming region is not exposed to an EUV/soft X-ray flux expected from an 
optically thick boundary layer. Kallman's (1983) photoionization analysis of 
the abundances of Si IV, CIV, and NV in the ultraviolet line forming region of 
the dwarf nova TW Virginis in outburst also results in the conclusion that a 
strong EUV/soft X-ray boundary layer flux cannot be illuminating the gas 
responsible for the ultraviolet lines. 
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These results should be regarded with caution, however, Ferland et al. 
assume a constant density n e = 10 10 cm" 3 . This Is probably a poor 
assumption. For example, King et_jrt_. (1983) find n e > 10* 2 cm" 3 for the line 
forming region In the nova-like variable UX UMa. It Is probable that there 
are large density gradients In the clrcumstellar gas of CVs, with n e > 10 12 
cm" 3 at the base of a wind and n e much smaller farther out In the wind. 

Kallman assumes that the Si IV, C IV, and N V absorption features are all 
formed in the same region. If this is not the case, his conclusions must be 
modified. For UX UMa, King et al_. (1983) find that the C IV, Si IV, and N V 
lines have different eclipse behavior, suggesting that the lines are not 
formed In Identical regions. As K1ng_et^K note, a fully self-consistent 
model for the clrcumstellar gas In CVs Is required to constrain the nature of 
the boundary layer radiation. 

Ferland jat_al_. (1982a) discussed the possibility that X-ray emission is 
suppressed by absorbing gas in nova systems. They concluded that there is not 
sufficient absorbing gas in the corona of V603 Aquilae to suppress the X-ray 
emission and that this scenario was therefore an unattractive explanation for 
the "nova problem". This conclusion may not be correct. While the column 
densities of 10 2 * - 10 22 cm" 2 inferred for the V603 Aquilae model corona of 
Ferland _et_al_. (1982b) are not sufficient to absorb a 1 - 10 keV hard X-ray 
flux, they are sufficient to cut off a 0.1 - 0.5 keV soft X-ray flux unless 
most of the helium is doubly ionized. 

McCray and Lamb (1976) have proposed the existence of an opaque shell 
around the neutron star in Her X-l to explain the strong soft X-ray flux 
observed from Her X-l as reprocessed flux from the central hard X-ray 
source. If a similar process occurs in CV S , in which a central soft X-ray 
flux irradiates an optically thick region which reprocesses it into longer 


wavelengths, boundary layer models can be reconciled with X-ray 
observations. If, on the other hand, It can be demonstrated that such a 
process Is unlikely, a reconsideration of boundary layer theory may be 
required. 

c) Other Types of CVs 

The preceding analysis of HEAO-1 constraints on the X-ray emission from 
novae cannot be used for other CV subclasses. The nova subclass Is uniquely 
suitable for an analysis constraining the fraction of the accretion luminosity 
radiated by an X-ray component because the number density of the novae Is a 
function of the white dwarf mass and mass accretion rate, through Its 
dependence on the outburst recurrence time scale (§IIIb). 

The Einstein IPC observations of selected novae, nova-like variables, and 
dwarf novae may help us to understand the nature of their boundary layers, 
when considered in the context of the constraints on the boundary layers in 
novae. Particularly useful is the ratio of the received flux from the hard 
X-ray component in the 0. 1-4.0 keV band, f x , to the received flux in the 
visual band, f v . Mason and Cordova (1982) note that the dwarf novae in 
quiescence have larger f x /f v ratios than the dwarf novae in outburst. 

Although there is more scatter in f x /f v for the novae and nova-like variables, 
as a group they tend to have f x ./f v ratios closer to those of dwarf novae in 
outburst than dwarf novae in quiescence. 

It is more difficult to use f x /f v ratios to test boundary layer models 
than it is to use, ratios of the X-ray component luminosity to the total 
accretion luminosity, because the ratio of visual band flux to total 
bolometric flux may depend on poorly determined factors such as reprocessed 
fluxes from the white dwarf, the accretion disk, and the secondary star. It 


Is tempting though to infer from the similarity in f x /f v between the novae, 
the rova-Hke variables, and the dwarf novae In outburst that nova-like 
variables and dwarf novae in outburst have boundary layers similar to nova 
boundary layers. Since the dwarf novae in quiescence are believed to have 
lower mass accretion rates, it may be that they produce hard X-rays In an • 
optically thin boundary layer while the boundary layers In novae, nova-like 
variables, and dwarf novae In outburst are optically thick soft X-ray 
producers. 

Subclasses of CVs which are not believed to have accretion disk boundary 
layers, such as the AM Herculls stars, have not been considered in this work. 

d) The Hard X=Ray Source 

If novae, nova-like variables, and dwarf novae in outburst typically have 
optically thick boundary layers, what is the source of the v/eak hard X-ray 
flux observed from them? The HEAO-1 constraints cannot directly Identify this 
source, but they do constrain the fraction of the accretion luminosity which 
can be radiated by this source for the novae and, by Inference, for the 
nova-like variables and dwarf novae in outburst. For example. If we adopt N H 
- 10 2 * cm" 2 as typical, Figure 5 shows that less than }% of the accretion 
luminosity can be radiated by a 1-10 keV optically thin plasma. Even for 
as large as 10 22 cm“ 2 , less than 3 % can be radiated by the hard X-ray source. 

Where in a CV system can we find a hot thin plasma which radiates only a 
small fraction of the gravitational energy released by accretion? Jensen et 
al. (1983) have suggested that the weak hard X-ray flux observed from the 
nova-like variable TT Arietis is produced in an accretion disk corona. If 
only a small fraction of the energy dissipated in the accretion disk and 
boundary layer is transported to the corona, it can explain why only a small 
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fraction of the accretion luminosity Is radiated by the hard X-ray source. 

This explanation may apply to the weak hard X-ray fluxes generally observed 
from novae, nova-like variables, and dwarf novae In outburst. 

e) Concluding Remarks 

The lack of detection of novae In the HEAO-1 LED! survey, and the 
observation of very weak hard X-ray fluxes from novae, nova-1 Ike variables, 
and dwarf flovad In outburst by the Einstein IPC, can be consistent with 
boundary layer models for the accretion process in cataclysmic variables if 
novae, nova-like variables, and dwarf novae in outburst typically have 
optically thick boundary layers with boundary layer temperatures kT « 15-30 eV 
and neutral hydrogen column densities > 10 20,5 -10 21> cm" 2 . 

If so, the boundary layer Is the source of a strong soft X-ray flux which 
Is highly absorbed by circumstellar gas and re-radiated at longer 
wavelengths. A fully self-consistent model for the circumstellar gas In CYs 
Is desirable to Investigate this possibility. 

The weak hard X-ray fluxes observed from these systems are unlikely to be 
produced In the boundary layer, but may be produced in an accretion disk 
corona. 
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APPENDIX 

The number of nova outbursts N 0 which will ocuur during a time Interval 
t q at a distance < d 0 Is given by 


T. d o 2* 

N 0 * —■ / dr / rde 
r o o 


(d 2 - r 2 ) 1 ^ 2 

/ n(r,e,z) dz 


(Al) 


where t r Is the characteristic time scale for the recurrence of a nova 
outburst in a single nova system and n(r,e,z) Is the number density of the 
novae as a function of cylindrical coordinates centered at the sun. The 
ratio t q / tp Is the fraction of novae which are expected to have outbursts 
during a time interval t q . The Integration of the number density over a 
sphere with radius d Q centered at the sun gives the total number of novae 
within a distance d 0 

O O 1/2 

d 2 tt (d - r^) 

N(d ) « / 0 dr / rde / 0 n(r,e,z) dz. (A2) 

0 o*o o 


If we define d x as the maximum distance at which thf? X-ray emission from 
a nova can be detected by the HEAO-1 LEDl detector, the integration of the 
number density over a sphere with radius d x centered at the sun gives us the 
number of nova systems <N X > detectable by the HEAO-1 soft X-ray survey: 


<N X >= N(d x ) = / x dr 


2rr 

/ rde 
o 


2 2, 1 / 2 

(d -r 2 ) 

/ n(r,e,z) dz. 

o 


(A3) 


In the CJN and Becker treatments, n(r,e,z) - constant = 10"^ pc“ 3 . A more 
careful analysis should account for scale heights in the distribution of the 
novae. For my analysis, I use: 


T'vrjr 




•• - ~ **■*#*#& • • - r rwW‘ ■ 




26 


n(p,+,*> - A8- (p/3 - 4 k P c > o- (z/0 - 2 k P c) 


(A4) 


where (p,$,z) are cylindrical coordinates centered at the galactic center and 
A Is an arbitrary constant. The scale heights h rt * 3.4 kpc and h z * 0.2 kpc 
are estimated from the distribution of gas and dust In the galaxy (Allen 
1973). Transforming n'p,<j>,z) In coordinates centered at the galactic center 
to n(r,,6,z) In coordinates centered at the sun (p «10 kpc, z ■ 0), we derive 


n(r,Q,z) » Ae 


-[(100 + r 2 + 20 r ccse) 1/2 /3.4] -5z 

6 


(A5) 


where r and z are in units of kpc. 

Combining equations Al, A3, and A5, we derive: 


<N X > T r I(d x ) 


(A6) 


where: 


1(d) a /dr /de [r e-^M^Orcose) 1 ' 2 ^,, . e -5(d 2 -r 2 ) 1 / 2 ) (A7) 
0 0 


and d Is In units of kpc. See Jensen (1982) for details. 

Using the estimates N 0 - 12, t q s 88 years, d 0 = 1 kpc, we have: 


<V T r I(d x ) 
T2 = BTyears TTTTOT 


(A8) 


Since 1(1.0) = 0.1518, we derive: 


1.1 <N V > 

I( V = WZH year) 


(A9) 
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which Is the equation for d x used in §IIIb. 
From equations Al, A5, and A7, we have: 


No ■ T- A 1 (d o> 
l r 


(A10) 


If d Q Is in units of kpc, A is in units of kpc” 3 . From our estimates N 0 = 
12, t q * 88 years, and d 0 = 1 kpc, hence I(d 0 ) = 0.1518, we can estimate A: 


A ■ 9 * 10 ‘ 9 <tW 


pc 


-3 


(All) 


and n: 


n(p,<|>,z) * 9xl0“ 9 


fir ,\ e -p/3.4 kpc 

Vl year' 


0 -z/O.2 kpc 



(A12) 


The estimated number density of the novae is proportional to the outburst 
recurrence time scale because we estimate the number density of all novae from 
the number density of observed novae. We must therefore multiply the observed 
number density by the inverse of the fraction of all novae which have been 
observed. This inverse fraction is estimated as 

In the neighborhood of the sun (p « 10 kpc, z « 200 pc), the number 
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density of the novae Is then estimated to be 


n * 4.7 x !0“ 10 ( 


r year 


) pc 


-3 


(A13) 


This relation Is reflected In Figure 1 
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FIGURE CAPTIONS 


Figure 1 - The dependence of x r% the timescale for the recurrence of a nova 
outburst, on the white dwarf mass and mass accretion rate, according to 
the models of Narlal and Nomoto (1979). The local space density of 
novae, n, which Is proportional to t r , Is also Indicated. 

Figure 2 - The function 1(d) Is shown for distances between 10 pc and 2 
kpc. 1(d) is a measure of the relative number of novae located at 
distances < d from the sun. The estimated vertical scale height in the 
distribution of novae (h z = 200 pc) Is Indicated. For d « h z , I a d 3 ; 
for d » h z , lad 2 . 

Figure 3 - The dependence of d x , the maximum distance at which a nova Is 
detectable by the HEA0-1 soft X-ray survey, on the white dwarf mass (M) 
and mass accretion rate (M). Since It is derived under the assumption 
that one nova Is near enough to be detectable by HEA0-1, d x is also the 

p * 

expected distance to the nearest nova. The trend d x a M is Indicated by 

% 

the dashed line. 

Figure 4 - The dependence of $ s ( GMM/R)“^ 4ir d x 2 (4.8 counts s" M (10" 11 
ergs cm" 2 count" 1 ) on the white dwarf mass (M) and mass accretion 
rate (M). The relatively small variation in <f> over the wide range of M 
and M means that the fraction of the total accretion luminosity which can 
be radiated as an X-ray component by novae, f, can be well constrained 
independent of mass and mass accretion rate. 

Figure 5 - Upper limits to f, the fraction of the total accretion luminosity 
that can be radiated by novae as a thermal bremsstrahlung spectrum, as a 
function of the temperature (kT) and neutral hydrogen column density N H 
of the spectrum. Contours of constant N H are shown. The log^g (cm" 2 ) 
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Is Indicated for each contour. The horizontal dashed line marks 
f « 0.25. The HEAO-1 observations are consistent with a boundary layer 
source only for thermal bremsstrahlung spectra characterized by 
temperatures and column densities lying above the dashed line. 

Figure 6 - Upper limits to f r the fraction of the total accretion luminosity 
that can be radiated by novae as a blackbody spectrum, as a function of 
the temperature (kT) and neutral hydrogen column density (N H ) of the 
spectrum. Contours of constant are shown. The log^g (cm" 2 ) Is 
Indicated for each contour. The horizontal dashed line marks f < 0.25. 
The HEAO-1 observations are consistent with a boundary layer source only 
for blackbody spectra characterized by temperatures and column densities 
lying above the dashed line. 

Figure 7 - The fraction of the 180-280 eV energy flux emitted from a 30 eV 
blackbody which is received after passing through the indicated column 
densities. This fraction becomes very sensitive to column densities 
greater than ~ 10 20 cm” 2 . 
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